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We present results of full 3D hydrodynamical and radiative transfer simulations of the collid¬ 
ing stellar winds in the massive binary system 77 Carinae. We accomplish this by applying the 
SimplcX algorithm for 3D radiative transfer on an unstructured Voronoi-Delaunay grid to re¬ 
cent 3D smoothed particle hydrodynamics (SPH) simulations of the binary colliding winds. We 
use SimpleX to obtain detailed ionization fractions of hydrogen and helium, in 3D, at the reso¬ 
lution of the original SPH simulations. We investigate several computational domain sizes and 
Luminous Blue Variable primary star mass-loss rates. We furthermore present new methods 
of visualizing and interacting with output from complex 3D numerical simulations, including 
3D interactive graphics and 3D printing. While we initially focus on 77 Car, the methods em¬ 
ployed can be applied to numerous other colliding wind (WR 140, WR 137, WR 19) and dusty 
‘pinwheel’ (WR 104, WR 98a) binary systems. Coupled with 3D hydrodynamical simulations, 
SimpleX simulations have the potential to help determine the regions where various observed 
time-variable emission and absorption lines form in these unique objects. 


1 3D SPH Simulations of rj Car’s 
Colliding Winds 

As already discussed in these proceedings (e.g. Gull, 
Morris et ah, Corcoran et ah, Hamaguchi et ah), 
Eta Carinae (77 Car) and its surrounding bipolar Ho¬ 
munculus nebula constitute an ideal astrophysical 
laboratory for the study of massive binary interac¬ 
tions and evolution, and stellar wind-wind collisions. 
Recent three-dimensional (3D) simulations set the 
stage for understanding the complex time-dependent 
3D stellar wind outflows in 77 Car. Spurned by re¬ 
cent suggestions that the mass-loss rate of the Lumi¬ 
nous Blue Variable (LBV) primary star might have 
dropped by a factor of 2 — 3 over roughly the past 
decade (Corcoran et al. 2010; Mehner et al. 2012), 
Madura et al. (2013) presented results from a series 
of 3D smoothed particle hydrodynamics (SPH) sim¬ 
ulations of 77 Car’s binary colliding winds assuming 
three different primary star mass-loss rates and us¬ 
ing various computational domain sizes. 

These SPH simulations reveal that the value of 
the primary’s mass-loss rate greatly affects the time- 
dependent hydrodynamics of the wind-wind colli¬ 
sion at all spatial scales investigated (~15.5 AU to 
1550 AU). The SPH simulations also show that the 
post-shock wind of the companion star switches from 
the adiabatic to the radiative-cooling regime dur¬ 
ing periastron passage {(j) « 0.985 — 1.02). This 
switchover starts later and ends earlier the lower the 
value of the primary’s mass-loss rate and is caused by 
the encroachment of the primary wind into the ac¬ 
celeration zone of the companion’s wind, plus radia¬ 
tive inhibition of the companion’s wind by the more 
luminous primary. This cooling switchover has im¬ 


portant implications for understanding the peculiar 
behavior of 77 Car’s X-ray light curve during the past 
two periastron passages, which showed an inconsis¬ 
tent recovery from the X-ray minimum (Corcoran et 
al. 2010, Corcoran et al., these proceedings). 

Understanding the behavior of the post-shock 
companion wind in 77 Car is crucial for understanding 
the behavior of the X-ray light curve, since it is in the 
extremely hot (up to ~T0 8 K) post-shock secondary 
wind that the majority of the hard X-rays observed 
are generated. As discussed in detail in Madura et 
al. (2013), the varying duration of 77 Car’s X-ray min¬ 
imum from cycle-to-cycle likely depends strongly on 
the length of time the post-shock secondary wind 
can remain in the radiative cooling regime during 
periastron passage. For most of 77 Car’s highly ec¬ 
centric orbit (e ~ 0.9), the post-shock companion 
wind cools adiabatically and thus remains relatively 
hot, generating hard X-rays. However, the cooling 
parameter \ that determines whether the post-shock 
gas is in the adiabatic or radiative cooling regime 
is proportional to the pre-shock wind speed to the 
fourth power (Stevens et al. 1992). Thus, even a 
small change in the pre-shock wind speed of the com¬ 
panion star can have a dramatic effect on the cooling 
of the post-shock gas. 

Normally, we think of the radiatively-driven wind 
from a massive star as following a so-called /3-law, 
wherein the wind speed at radius r is given by 
v(rj = 77oo(l — R*/r)P , where Voo is the wind termi¬ 
nal speed and i?* the stellar radius. However, as dis¬ 
cussed in Madura et al. (2013), because of the highly 
elliptical orbit and presence of the incredibly lumi¬ 
nous primary star, the pre-shock companion wind 
speed in the 77 Car system is a complicated function 
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of the stellar luminosity ratio and stellar separation. 
The pre-shock wind speed of the companion star in 
rj Car decreases by roughly a factor of three between 
apastron (phase 0.5) and phase 0.99 (Madura et al. 
2013). Because the cooling parameter y depends on 
the wind speed to the fourth power, a factor of three 
reduction in the companion’s wind speed around pe- 
riastron causes the post-shock gas to switch from 
being adiabatic to strongly radiative. This reduces 
greatly the volume of hot X-ray generating gas, and 
thus the observed X-ray flux. Only after periastron 
when the stellar separation starts to increase does 
the pre-shock companion wind speed also increase, 
and the post-shock gas can switch back to the adi¬ 
abatic regime, again producing hard X-rays. The 
duration of the radiative cooling phase of the post¬ 
shock companion wind thus very likely affects the 
duration of the observed X-ray minimum and tim¬ 
ing of the X-ray recovery following periastron. 

Since the primary star in r) Car is a known LBV 
subject to periods of instability, and because the 
wind-wind collision region is very chaotic and dy¬ 
namic due to the various instabilities that arise 
at the interface between the two colliding winds 
(Stevens et al. 1992; Madura et al. 2013), it is reason¬ 
able to suggest that such instabilities and/or small 
variations in stcllar/wind parameters can lead to 
changes in the overall pre-shock wind speeds, loca¬ 
tion of wind momentum balance, and duration of 
strong radiative cooling in the post-shock compan¬ 
ion wind. Thus, changes in the duration of the X-ray 
minimum and timing of the X-ray recovery following 
periastron may be natural in a system like rj Car. 

2 3D Radiative Transfer Models 

In addition to the ‘current’ colliding winds interac¬ 
tion responsible for the X-ray emission, which occurs 
in the innermost regions of the rj Car system (at 
spatial scales comparable to the orbital semimajor 
axis length a = 15.4 AU), larger scale (~ 3250 AU 
diameter) 3D SPH simulations reveal outer wind- 
wind collision regions (WWCRs) that extend thou¬ 
sands of AU from the central stars (Madura et al. 
2012, 2013). Long-slit spectral observations with 
the Hubble Space Telescope/Space Telescope Imag¬ 
ing Spectrograph (HST/STIS) reveal these spatially- 
extended WWCRs, seen via emission from multiple 
low- and high-ionization forbidden lines (Gull et al. 
2009, 2011; Teodoro et al. 2013). Observations of 
the different broad high- and low-ionization forbid¬ 
den emission lines provide an excellent tool to con¬ 
strain the orientation of the system, the primary’s 
mass-loss rate, and the ionizing flux of the hot sec¬ 
ondary star. In addition to emission from forbid¬ 
den lines, helium presents several interesting time- 
variable emission and absorption features that pro¬ 
vide important clues on the geometry and physical 
proprieties of the system and the individual stars. 


To set the stage for future efforts aimed at generat¬ 
ing synthetic observations for comparison with avail¬ 
able and future HST/STIS and ground-based data, 
we have performed a series of 3D radiative trans¬ 
fer simulations of the interacting winds in 77 Car. 
We use the SimpleX algorithm to post-process the 
3D SPH simulation output and obtain detailed 3D 
maps of the ionization fractions of hydrogen and he¬ 
lium. The resultant ionization maps constrain the 
regions where observed forbidden emission lines can 
form, and the regions where helium is singly- and 
doubly-ionized. 

Our SimpleX results recently appeared in the se¬ 
ries of papers Clementel et al. (2014, 2015a,b). In 
Clementel et al. (2014), we focus on large-domain 
(~3250 AU diameter) simulations at an orbital 
phase of apastron, investigating the 3D ionization 
structure of the extended colliding wind structures 
assuming different primary mass-loss rates. These 
simulations are comparable in size to the HST/STIS 
observations of the extended forbidden emission 
structures (Gull et al. 2009, 2011; Teodoro et al. 
2013). We find that at apastron, the dense primary 
wind confines the photoionization, and hence the for¬ 
bidden line emission, to the secondary star’s (and ob¬ 
server’s) side of the system. Changing the primary’s 
mass-loss rate results in quite different ionization 
volumes, with much more ionized gas present for 
lower mass-loss rates. The large apparent changes in 
ionization volume with decreasing primary mass-loss 
rate imply that any major decrease in primary mass- 
loss should lead to significant observable changes in 
the spatial extent, location, and flux (proportional 
to density squared) of the broad high-ionization for¬ 
bidden emission lines. Future models based on the 
SimpleX results may be used to constrain any such 
potential mass-loss changes. 

Clementel et al. (2015a) and Clementel et al. 
(2015b) investigate the ionization structure of he¬ 
lium in 77 Car’s innermost regions (<310 AU in di¬ 
ameter) at apastron and periastron, respectively, for 
different primary mass-loss rates. We find that re¬ 
ducing the primary’s mass-loss rate increases the vol¬ 
ume of He + . Doubly-ionized helium (He ++ ) only ex¬ 
ists in the extremely hot post-shock secondary wind. 
At orbital phases near apastron, lowering the pri¬ 
mary mass-loss rate produces large variations in the 
volume of He + in the pre-shock primary wind lo¬ 
cated on the periastron side of the system. Our re¬ 
sults show that binary orientations in which the sec¬ 
ondary star at apastron is on our side of the system 
are most consistent with the available observations. 

At an orbital phase of periastron, our simulations 
show that generally, He 0+ -ionizing photons from the 
companion star are not able to penetrate into the 
pre-shock primary wind, and that He + due to the 
companion’s ionizing flux is only present in a thin 
layer along the leading arm of the inner WWCR 
and in a small region close to the stars. Lowering 
the primary’s mass-loss rate allows the companion’s 
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ionizing photons at periastron to reach the expand¬ 
ing unshocked secondary wind on the apastron side 
of the system, and create a low fraction of He + in 
the pre-shock primary wind. With the companion 
at apastron on our side of the system, our results 
are qualitatively consistent with the observed varia¬ 
tions in strength and radial velocity of 77 Car’s he¬ 
lium emission and absorption lines. 

3 3D Visualization and Printing 

3D hydrodynamical and radiative transfer simula¬ 
tions of 77 Car show that the complex time-varying 
WWCR has a major impact on the observed X-ray 
emission, the optical and ultraviolet light curves and 
spectra, and the interpretation of various line pro¬ 
files and interferometric observables (see Madura et 
al. 2015, and references therein). While 3D simula¬ 
tions have helped to increase substantially our un¬ 
derstanding of the present-day 77 Car binary, histor¬ 
ically we have been limited by an inability to ade¬ 
quately visualize the full 3D time-dependent geom¬ 
etry of the WWCR. However, this is crucial if we 
are to thoroughly understand how and where various 
forms of observed emission and absorption originate. 

In an effort to better understand the 3D struc¬ 
ture of 77 Car’s time-dependent WWCRs, Madura 
et al. (2015) recently published 3D interactive visu¬ 
alizations and 3D prints of 77 Car’s colliding stel¬ 
lar winds at orbital phases of apastron, perias¬ 
tron, and three months after periastron. The fig¬ 
ures in the PDF version of Madura et al. (2015), 
downloadable directly from the NASA ADS service 
(http://adsabs.harvard.edu/), are fully 3D in¬ 
teractive when using the standard Adobe Reader 
(https://get.adobe.com/reader/). The 3D mod¬ 
els presented in Madura et al. (2015) also constitute 
the first 3D printable output from a supercomputer 
simulation of a complex astrophysical system. 

The 3D visualizations and prints in Madura et al. 
(2015) reveal important, previously unknown ‘finger¬ 
like’ structures at orbital phases shortly after perias¬ 
tron (~ 1.045) that protrude radially outwards from 
the spiral WWCR. We speculate that these fingers 
are related to instabilities (e.g. thin-shell, Rayleigh- 
Taylor) that arise at the interface between the radia- 
tively cooled layer of dense post-shock primary-star 
wind and the faster (~ 3000 km/s), adiabatic post¬ 
shock companion-star wind. The 3D visualizations 
also reveal that the inability of the companion star’s 
wind to collide with the prmiary’s downstream wind 


produces a large ‘hole’ in the trailing arm near the 
WWCR apex at periastron. The size of this hole 
is directly connected with the wind momentum ra¬ 
tio, and therefore the WWCR opening angle. As ex¬ 
pected, higher primary mass-loss rates cause a larger 
hole than lower primary mass-loss rates. After pe¬ 
riastron, this hole expands and propagates down¬ 
stream along the vanishing WWCR trailing arm. 

The success of the work in Madura et al. (2015) 
and the relatively easy identification of previously 
unrecognized physical features highlight the impor¬ 
tant role 3D printing and interactive 3D graph¬ 
ics can play in the visualization and understand¬ 
ing of complex 3D time-dependent numerical sim¬ 
ulations of astrophysical phenomena. The 3D print¬ 
able models of 77 Car’s colliding winds and the 77 Car 
Homunculus nebula (Steffen et al. 2014) are freely 
available online at the NASA 3D Resources website 
http://nasa3d.arc.nasa.gov/. 
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